Microscopic structure and dynamics of high and low density 
trans- 1 , 2-dichloroethylene hquids 
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We present a study of the dynamics and structural changes for trans- 1, 2-dichloroethylene between 
high and low density liquids using neutron scattering techniques (diffraction, small angle neutron 
scattering and time of flight spectroscopy) and molecular dynamics simulations. We show that 
changes in the short range ordering of molecules goes along with a change of the molecular dynamics: 
both structure and dynamics of the high density liquid are more cooperative than those of the low 
density liquid. The microscopic mechanism underlying the cooperative motions in the high density 
liquid has been found to be related to the backscattering of molecules due to a strong correlation of 
molecular ordering. 



Classical thermodynamics establishes the existence of 
one unique liquid state for any material, i.e., there is 
only a liquid phases characterized by its density given the 
thermodybamic coordinates pressure and temperature. 
Nevertheless, recent experimental results and molecu- 
lar dynamics (MD) simulations suggest that, even for 
one-component systems, several liquid phases can ap- 
pear with an associated liquid-liquid phase transition 
(LLPT). A noticeable number of cases has been found 
for atomic liquids, the best-known example concerning 
liquid phosphorus [l|, 0|, where the LLPT appears as 
a transition between thermodynamically stable phases 
with strong structural changes 3j]. As far as molecular 
liquids are concerned, the number of experimental evi- 
dences for LLPT is still rather scarce and comprises only 
a limited number of compounds such as triphenyl phos- 
phite m and n-butanol [5|. According to the so called 
two order parameter theories that propose an explana- 
tion for the LLPT, liquids must be described not only 
by their density but also by an additional order parame- 
ter accounting for changes in the molecular arrangement 
[^-Q . The LLPT can end in a liquid- liquid critical point 
between a high-density liquid (HDL) and a low-density 
liquid (LDL). However, changes in the dynamics with 
an associated change in structural featurjjs can also be 
explained by a singularity- free scenario (^]. In the lat- 
ter case, changes in both dynamics and structure from a 
HDL to a LDL also take place at the point where the iso- 
baric heat capacity Cp has a maximum, but no critical 
point or LLPT are observed at non-zero temperature. 

An early work on trans-1, 2-dichloroethylene (Tmoit — 



223 K) suggested the existence of a LLPT at Tt = 
247 K > Tnicit based on a small jump in density (less than 
0.06%) as well as in the compressibility, and a clear dis- 
continuity on the spin-lattice relaxation time Ti l3, U | ■ 
The observed changes in Ti were tentatively related to a 
lack of freedom of the molecular rotation in the HDL, not 
present in the LDL. The change on the dynamics of this 
substance between both liquids, was thereafter also sup- 
ported by discontinuities in the viscosity measurements 
and the slope of the rotational relaxation time [12], and 
by the absorbance, frequency and linewidth of several 
infrared vibrational spectroscopy bands (13j . Concerning 
structural related magnitudes, subsequent measurements 
of the density as a function of the temperature did not 
find a jump at the expected LLPT but only a change in 
its slope [12] ■ More recently, some of the authors of the 
present work have also undertaken calorimetric and non- 
linear dielectric experiments [l4{. In that work a strong 
pre-transitional anomaly of nonlinear dielectric effect was 
obtained, similar to the one observed in the isotropic 
phase of nematic liquid crystals, together with a max- 
imum of Cp. Therefore, experimental results unambigu- 
ously show a clear change in the dynamics, together with 
a slight change in the structure, between HDL and LDL, 
that takes place when Cp is in a maximum. However, 
those facts are not enough to unambiguously determine 
if they are related to a singularity-free or a liquid-liquid 
critical point scenario |9|, i.e., the liquid undergoes a first 
order phase transition. The present study is aimed to in- 
vestigate the microscopic structural and dynamical differ- 
ences between the HDL and LDL, from the experimental 
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FIG. 1: (color online). Total interference function of trans- 
1,2-dichloroethylene for the LDL and the HDL, obtained by 
neutron scattering experiments (a) and by molecular dynam- 
ics simulations (b), where an ab initio calculation of the 
molecular structure is also shown. Inset shows the macro- 
scopic cross section for low Q values obtained by SANS ex- 
periments also for temperatures above and below Tt. 



point of view and from MD simulations. 

Because the differences between the HDL and LDL in 
trans- 1,2-dichloroethylene were related to a change from 
a nematic-like to an isotropic molecular ordering as tem- 
perature is raised [13, Ell, we have performed a series 
of small angle neutron scattering (SANS) measurements 
from 220 K to room temperature to ascertain whether 
there is formation of intermediate range nematic-like 
structures or clustering in the HDL. Experiments were 
performed on the KWS-2 diffractometer of the Jiilich 
Centre for Neutron Science at the Forschungsneutronen- 
quelle Heinz Maier-Leibnitz (FRM II, Munich, Germany) 
[15i] using a wavelength of 4.5 A and a sample-detector 
distance of 2.0 m that allowed to perform measurements 
in the Q range between 0.003 and 0.3 A-\ The pro- 
gram QTIKWS iJJj] was used to perform data correction 
and normalization. Results are shown in the inset of Fig. 
[l]for the HDL and LDL. No pronounced signal has been 
obtained within the experimental error in the measured 
Q-range, which disfavors the existence of a long range 
nematic-like ordering for the HDL The data exclude 
that the differences between both liquids are related to 
the emergence of molecular clustering on length scales of 
about 20-2000 A. 

Going down in the spatial range, the microscopic short 



range order (SRO) concerning a length scale of the or- 
der of a few molecular lengths (/moi ~ 4 A) has been 
analyzed by means of neutron scattering experiments on 
the D4c diffractometer at the Institute Laue-Langevin 
(ILL, Grenoble, France) [3| using a wavelength of 0.5 A 
and a deuterated sample. Data were corrected and nor- 
malized using the program CORRECT [l^ and inelastic 
corrections were also carried out (for details on data re- 
duction see Re f. 12011 ). The obtained total interference 
function F{Q) [2l\ is shown in Fig. [T] for two representa- 
tive temperatures. A change in the shape of the profile 
emerges between the first and second peak which reflects 
a change in the SRO . This change was also observed in a 
series of temperature dependent experiments on the D20 
diffractometer (ILL) using a wavelength of 1.3 A, giving 
better access to the low-Q region. It should be pointed 
out that similar changes have been found in the case of 
experiments performed on HDL water at hig h pressures 



23|,l24|. 



[22| and MD simulations on HDL silicon 

The microscopic mechanisms giving rise to the changes 
in the interference function F{Q) have been investigated 
through a series of MD simulations (Fig, [p. Those were 
performed using the program AMBErS [25|,|26[ with a sim- 
ulation box containing 3629 molecules and a time step of 
1 fs. The inter- and intra-molccular interactions for the 
trans- 1 ,2-dichloroethylene molecule were described using 
the GAFF force field 



26l |. The simulations were car- 
ried out for the temperature range 200-350 K in the NPT 
ensemble, therefore allowing the box size to change, and 
the total time scale of each simulation run was 40-50 ns 
27 , 2^ . As can be seen in Fig. [U the agreement between 
simulations and experiment is excellent. Note that not 
only the shapes of the simulated F{Q) closely resemble 
the experimental ones, but also positions and intensities 
are equal within the experimental error. 

To emphasize quantitatively which are the SRO 
changes between the HDL and LDL, we show in Fig. 
[2] the probability of finding two molecules with a deter- 
mined relative orientation for increasing distances. To 
present the distance dependent SRO, the molecular co- 
ordination number (MCN) has been chosen rather than 
distance in order to avoid trivial effects due to density 
changes. Figure [5] shows that the relative orientation of 
nearest neighbors is virtually the same and strongly de- 
fined for both liquids, in such a way that their C-C vec- 
tors are parallel, i.e., with cos a = ±1. However, upon 
increasing distance between the molecules (third to tenth 
neighbor), a clear difference between the LDL and HDL 
can be seen in the SRO. For instance, for MCN between 
eight and ten, molecules are arranged in a orthogonal 
way for HDL and are randomly oriented for the LDL. 
Further analysis (not shown) tell us that besides the rel- 
ative orientation of the C-C axes of two molecules, the 
reported changes are also refiected in the relative posi- 
tion of two molecules and in the relative orientation of 
the planes defined by the Cl-C-H bonds of two molecules. 
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FIG. 2: (color online). Short range order for the HDL (left) 
and the LDL (right), at 200 and 300 K respectively. The y- 
axis is the cosine of the angle a formed between the vectors 
defined by the C-C double bonds of two neighbor molecules 
and the x-axis is the number of molecules surrounding a cen- 
tral one. The 2;-axis represents the probability of finding two 
molecules at a distance determined by the MCN with a rela- 
tive orientation defined by the a angle. 



As proposed in the frame of two order parameter theories 
p-lsl, this liquid can therefore not be simply character- 
ized by its density, but also a parameter reflecting the 
SRO should account for the changes occurring between 
the two liquids. 

To ascertain the influence of the aforementioned struc- 
tural changes in the dynamics of the system, we have per- 
formed a series of quasielastic neutron scattering (QENS) 
experiments for temperatures ranging from 220 to 300 K, 
conducted with the TOFTOF spectrometer at FRM II on 
a hydrogenated sample. Spectra were measured using an 
energy resolution of 60 /leV and a wavelength of 6 A, and 
data reduction was performed using the program frida 
[29']. Two representative spectra at Q = 0.3 at the 
LDL and HDL are shown in Fig.[3I^a). 

To have a first insight of possible changes in the dy- 
namics, a stretched exponential exp(— i/r)*^ was fitted 
to the intermediate scattering function, obtaining a de- 
crease of the exponent /3 in the HDL, which is related 
to a broadening of the relaxation time distribution, i.e., 
an increase of the cooperativity of the molecular motion. 
A deeper analysis was performed through a careful fit of 
data to several models u sing a Bayesian approach with 
the program fabada [30|-|32|. In order to keep the num- 
ber of physical parameters describing the data to a min- 
imum, we performed the fits to the whole scattering law 
S{Q,uj). Model selection was performed calculating the 
maximum of the likelihood C'^^^ for each model. 

The first model used to describe the data is composed 
by a diffusion motion plus an isotropic rotation of the 
molecule [33[. In this way, the only physical parameters 
to describe the whole experimental scattering function 
were the translational and rotational diffusion coefficients 
as well as the radius of rotation R. This simple model is 
able to describe LDL data giving rise to a good quality fit 
and a radius of rotation almost independent of temper- 



ature {R « 1.72 A, see inset in Fig. |3l[a)), in agreement 
with the aforementioned NMR measurements [l^, 11 1. 
For the HDL an inability of this model to describe ex- 
perimental data (on quantitative grounds, a decrease on 
the C"^^^ of the fit), makes the isotropic model for ro- 
tation not valid to describe molecular rotation. This is 
also reflected on the sudden drop of R below Tj (inset in 
Fig. El^a)). Even when assuming a free diffusion model 
for rotations 3J1 and assuming an anisotropic rotation of 
molecules (S^, |3a] , models could not account for spectra 
obtained in the HDL. Only adding a confined motion to 
the previous model, data could be successfully described. 
For the latter model, S{Q,uj) can be expressed as: 

S{Q,Lo) = [A{Q)S{co) + {1 - A{Q)) ■ L{lo)] 

(E) Srot{Q,u})(SSdm{Q,uj)(E)R{Q,u}) 

where A{Q) is the elastic incoherent structure factor of 
the confined motion, L(uj) is a Lorentzian function ac- 
counting for a confined motion, SrotiQ,i^) is the com- 
ponent accounting for the molecular rotation, S'diff(Q, w) 
the component for molecular diffusion [33| . and R{Q,w) 
is the instrumental resolution. Agreement with HDL 
spectra is shown in Fig. El^a). To estimate the length 
scale at which the confined motion is taking place, A{Q) 
values have been fitted with a model of diffussion inside 
a sphere, yielding R — 1.91 ± 0.07 A and, more realisti- 
cally, a threedimensional Brownian oscillator 37[, yield- 
ing a mean squared displacement = 2.12±0.07 A. 
Moreover, the obtained A{Q) are almost independent of 
temperature for the HDL, which means that the length 
scale of the confined motion is roughly temperature in- 
dependent for this liquid. 

As previously performed for the microscopic structure 
investigation, the dynamics obtained from MD simula- 
tions have been analyzed as well. Agreement between 
translational diffusion activation energy determined us- 
ing the neutron experiments (Ea = 76 ± 1 meV) and MD 
simulations {Ea — 75 ±2 meV) confirms that the sim- 
ulation is indeed describing the dynamics of our system, 
and that we are actually simulating the liquid phase in 
all temperature range (see also inset of figure 3b) . The 
normalized velocity autocorrelation function (VACF) is 
shown in Fig. ^h). The fastest decay for the HDL is 
higher that that of LDL indicating that interaction with 
neighboring molecules takes before in that phase. As can 
be seen, there is a clear change in the dynamics of the 
two liquids, the VACF reaching negative values for HDL. 
Although interactions between neighbor molecules at the 
LDL seem not to impede molecular diffusion, for the HDL 
a well defined backscattering effect emerges, displaying a 
VACF with two minima characteristic to that of hydro- 
gen bonded systems Q . Such a change in the dyriamics 
is not expected to be due to temperature effects [39| . Ad- 
ditionally, a change of density of the system would not 
be able to explain this change in the dynamics since a 
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FIG. 3: (color online), (a) Fits of QENS spectra for Q = 
0.3 at 300 K (shifted upwards) to a model considering 
diffusion and rotation motions, and to a model with an ad- 
ditional confined motion for 220 K (solid lines). Dotted lines 
show the experimental resolution. Inset shows radius of ro- 
tation in function of temperature for the model with only 
diffusion and rotation (dashed lines are a guide to the eye), 
(b) VACF obtained by MD simulations as a function of time, 
for temperatures from 200 to 320 K in steps of 20 K (solid 
lines for the LDL and dashed lines for the HDL). The inset 
shows the obtained mean square displacement from simula- 
tions (solid lines for the LDL and dashed lines for the HDL). 



change as large as about 15% in density is needed to pro- 
duce an effect on the liquid dynamics in other systems 
[3^. Therefore, only a change in the SRO is able to ac- 
count for the change in the dynamics seen by experiments 
and simulations. 

The need for an additional confined motion at low tem- 
peratures to describe the QENS data goes along with a 
backscattering effect observed in the simulations. We 
can therefore assert that a cooperative molecular mo- 
tion is present for the HDL but not for the LDL, which 
agrees with other experimental results [l3| . The complex 
dynamics of the HDL agrees with a stronger molecular 
ordering present in this liquid, where on average there are 
many orthogonally oriented molecular pairs. This strong 
correlation is partially lost in the LDL, where molecular 
movements are due to non-coupled movements of diffu- 
sion and rotation. 

We have shown that the microscopic ordering of 
molecules and molecular dynamics are different between 
the HDL and the LDL. The changes in the dynamics 
are not due to temperature or density effects, but due 
to changes in the SRO: for the HDL there are molecules 
perpendicularly oriented that are randomly oriented in 
the LDL. This accounts for the change of molecular dy- 
namics from simple non-cooperative motions in the LDL 
to cooperative motions for the HDL. However, if these ef- 
fects are to be explained in the frame of a singularity-free 
scenario or a liquid-liquid critical point scenario associ- 
ated to a first order phase transition remains an open 
question. 

The authors would like to thank C. A. Angell, G. 
Franzese, G. Sese and T. Strassle for helpful discussions 
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